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'"^ . Abstract. Renormahzation group method is one of the most powerful tool to obtain 

I I approximate solutions to differential equations. We apply the renormalization group 

^ ' method to Hamiltonian systems whose integrable parts linearly depend on action 

f^ , variables. We show that the renormalization group method gives the same approximate 

j^». ' solutions as canonical perturbation theory up to the second order of a small parameter 

f\j I with action-angle coordinates. 
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r. ' 1- Introduction 

>'■ 

S^ . Dynamical systems written by differential equations are useful to understand temporal 

j^ ■ evolutions of the nature. Exact solutions to the equations are not always obtained 

because of non-integrability of systems, and naive perturbation often yields secular 
terms which prevent us from getting approximate but global solutions. Singular 
perturbation techniques |jl|], eg averaging methods, multi scale methods, matched 
asymptotic expansions and WKB methods, are available to construct global solutions. 
However, they provide no systematic procedures for general systems because we must 
select a suitable assumption about the structure of a perturbation series. 

Recently, renormalization group method is proposed 0, Q as a tool for global 
asymptotic analysis of the solutions to differential equations. It unifies the techniques 
listed above, and can treat many systems irrespective of their features. We apply the 
renormalization group method to Hamiltonian systems, and compare it with canonical 
perturbation theory 0, Q, which is one of the most developed perturbation theory for 
Hamiltonian systems. In this letter, we show that the renormalization group method 
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also unifies the canonical perturbation theory. That is, the former and the latter give 
the same solutions to equations of motion up to the second order of a small parameter. 
We use action-angle coordinates as they are suitable for perturbed Hamiltonian 
systems, and Hamiltonians are 

H{I,0) = Ho{I) + eHi{I,e), (1) 

where both / and are A^-dimensional vectors, the integrable part Hq is 

Ho{I) = u; ■ /, (2) 

and Hi{I, 6) is periodic with respect to each element of 0. 

We derive an approximate solution with naive perturbation in section |^, and then we 
renormalize secular terms to constants of integration in section ^. Finally, in section ^ we 
compare the renormalized solutions with solutions obtained by canonical perturbation 
theory. 

2. Naive Solution 



The equation of motion for the system (|1|) is 

da; 

d^ 

where x = {I, 6) is a 2A^-dimensional vector and the symbol {■ , -Ja, is Poisson bracket 
with respect to the subscript, in this case, x. We expand cc as a series of powers of e. 



— = {x, Hoix) + eH^{x)}^, (3) 



X = X 



(0) I .^(1) ^ .2^(2) 



ex^'^ + e'x^'^ + . . . , (4) 



and then equation (^ becomes 



^ (a;(0) + ex(i) + e^x^^^ + ..) = U,H, + eH, + e^^ • x^'^ | (a.^) + . . . , (5) 

where, in the right-hand-side, we substitute x^^"^ to x after the Poisson bracket has been 
operated. 

The solution to 0{e^) is 

/(o)=ao, 0(°)=a;t + /3o, (6) 

where A^-dimensional vectors olq and /3g are constants of integration. 
The equation of motion for O(e^) is 

^ = {^.H,}M% (7) 

and hence the solution to O(e^) is 

^^'^ = {X, S,{x)}^+ t {x, (//i(x))J^, (8) 

where we introduced the symbols 

X=(a,/3), OL = ao, f3 = uj t + /3q, (9) 
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{■)^ represents average over t, and 

S,{x)^ldt{H,{x)-{Hi{x))t)- 
The following relation was also used 

which is satisfied by arbitrary functions / and g that are periodic for 6 and (3. 
The equation to O(e^) is 

da;(2) 



(10) 
(11) 



dt 



{x , {H^ix) , S,{x)}^}^+{{x , S,{x)}^ , H,{x)}^ 
+ t {{X. H,{x)}^ , {H,{x))t] , 



(12) 



and the solution to O(e^) is 



x^'^ = {X , S2}^ + 1{{X, S^}^ , S,]^ + t {x, (F,),}^ 

+ (t^-secular terms) + (t-secular terms with non-constants). 



(13) 



Here 



and 



S2{x)^Jdt{F,{x)-{F2{x)),), 

F2{x)^{Hi{x) , S,{x)}^ + ^{{Ho{x) , S,ix)}^ , S,ix)} ^, 

{{X, S^}^ , hA =|x, ^{{^0, S,}^ , SA I +A{x) + B{x), 

JdtA{x) = l{{x, S^}^ , 5i}^, 

jdtt {{x, H^}^ , {H,),]^ = - jdtB{x) 



(14) 



+ (t^-secular terms) + (t-secular terms with non-constants), 

which are proven by using Fourier expressions of i^i(x), {Hi{x))t ^^^ Si{x)- 
concrete forms of the Fourier expressions, A{x) and B{x) are shown in appendix. 
Consequently, the naive solution to equation (H) is, up to O(e^), 



(15) 
(16) 

(17) 
The 



x = X + e{x , Si} +e'^ 



{X, ^2}^ + ^{{x, S,}^ , S,]^ 



+ t 



e{x, {H,),}^ + e'{x, (^2),}^ 



:i8) 



+ e^ [(t^-secular terms) + (t-secular terms with non-constants)] . 
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3. Renormalization of Secular Terms 



We renormalize the secular terms of the naive solution ([T8|) to the constants of 
integration. First we regard ctQ and /3q as functions of t which are 



(Xo{t)=<Xo + t e{a, {H^)^}^ + e^ {a , (Fs),! 



(19) 



/3o(t) =f3, + t[e{f3, {H,),}^ + e' {f3 , (F^) J^j 

Second we introduce assumptions with which the renormalization transformations (p!9|) 
becomes to be a Lie group. In this case, we assume that equation (p^Of ) is a truncated 
Taylor series of ao(t) and /3o{t) around the initial time t = . From time-evolutional 
symmetry of the system (Q), the renormalization group equation becomes 

dcKo 



dt 

dt 



O(e') 



in other words, 



dx 



dt '^ ' "^^-^^x 
The renormalized solution is therefore 



e{«, (/Ji)J^ + e^{«, (F2)J^ 
{X, ^o(x)L + e{x, (i/iU +eMx, (^2),L + 0(e 



(20) 



tJx 



'<Jx 



a; 



X + e{X , '^il^ + e' 



{X, 52}^ + ^{{X, ^ilx ' '^i}. 



+ 0(e' 



(21) 



(22) 



where x is governed by equation (plf ). Here, t^-secular terms and t-secular terms with 
non-constants in O(e^) of equation (p!8| ) are renormalized to coefficients oft-secular terms 
and coefficients of non-constant terms, respectively. 

4. Comparison with Canonical Perturbation Theory 

Finally we compare the renormalized solution (^2]) and solution obtained by canonical 
perturbation theory. The strategy of the theory is to canonically transform coordinates 
X = {I, 0) to X* = (/*, 0*) with the generator S{x*) 

X = exp{eDs)x*, Dsf{x*) = {f{x*) , S{x*)}^. , (23) 

such that secular terms do not appear in the coordinates x*. What we must calculate 
are the generator S and the transformed Hamiltonian H* . Canonical perturbation 
theory [^, |[ states that the required generator S = Si + eS'2 + ... is expressed as 



Si{x*) = I dr {Hi{x*) - {Hi{x*))^), 



(24) 
S2ix*) = j dr {F,{x*) - {F,{x*))^), (25) 

and the transformed Hamiltonian H* = Hq + eHl + e^i/g + ■ ■ ■ ^s 

H;{x*) = H,{x^), Hlix*) = {Hiix*))^, H;ix*) = {F,{x*))^, (26) 
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where 

F,{x*) = {H^{x*) , S,{x*)}^, + 1 {{Hoix*) , S,{x*)}^. , S,{x*)}^, . (27) 

The symbol (•)^ represents the average over r, the time of x* following H^, that is, 

dx* 

^ = {^* , ^o(^*)}.* • (28) 

Consequently, this theory gives an approximate solution determined by 
da;* 

and the canonical transformation 



{x* , Ho{x*)}^,+e{x* , {H,)^}^, + e'{x* , (F,) ^} ^, + O (e^) (29) 



X = x* + e{x* , Si}^, + e 



2 

X* 



K, S,}^, + ^-{{x\ S^}^. , S,}^, 



+ O(e=^).(30) 



The approximate solution ( pOD is the same as the renormalized solution (^) since 
temporal evolutions of x* and x ^i^re governed by equations (|29| ) and (|2T| ) respectively, 
and the two equations have the same structure. 

5. Summary and Discussions 

We showed that renormalization group method gives the same approximate solutions 
as canonical perturbation theory up to the second order of a small parameter to the 
Hamiltonian systems whose integrable parts linearly depend on action variables. That is, 
renormalization group method unifies not only averaging methods, multi scale methods, 
matched asymptotic expansions and WKB methods, but canonical perturbation theory. 
We suppose that the unification holds even in higher orders of the small parameter. 

In systems whose integrable parts are not linear, secular terms are not always 
proportional to time t, and may be proportional to t^ [n ^ 1). Canonical perturbation 
theory cannot remove the latter secular terms since subtracting time-averages of 
perturbative part of Hamiltonian is effective only for the t-linear secular terms. On 
the contrary, renormalization group method gives global solutions by introducing 
assumptions with which renormalization transformations become to a Lie group [^ 
and can treat t" type secular terms. 

In the previous paper [Q, we discussed relation between integrability of original 
systems and symplectic properties of renormalization group equations in Cartesian 



coordinates. From equation (|2T|) , we clarified that renormalization group equations 
are always Hamiltonian systems in action-angle coordinates whose Hamiltonian is 

H''''{x) = Ho + e{H^), + e'{F2),. (31) 

Symplectic properties are recovered even in Cartesian coordinates by using "gauge 
freedom" which is homogeneous terms of O(e^). Details will show in the next paper 0. 
YN is supported in part by the Grand-In-Aid for Scientific Research of the Ministry 
of Education, Science, Sports and Culture of Japan(09740196). 
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Appendix 

Let us introduce Fourier series of Hi, {Hi)^ and Si as 



^i(x) =Xl^i'^("^^ 



iml3 



(Hi), ix) = Yl HiAo.y-^-^ 

rri: Res 

Siix) = E -^—Hi 



m-.Rcs 

1 



im ■ uj 

TTiiNon 



where the symbols YlmKcs ^'^'^ XlrriNon represent to take summations over m such that 
m ■ a; = and m ■ a; 7^ 0, respectively. By using these expressions and r) = {—6, 1), 
the concrete forms of A{x) and B{x) are 






c{x) 



Tn,;Non riiNon 

and 



m ■ u) 

miRcs TiiNon 



•^— ' •^— ^ in . C.I 

respectively, where 
C{x 



. . dHi,^ , d{in -13) f^ \ dHi 



or) or) J oak 

i{m,+n)-f3 



dakdr) dr) dak 



e 



